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We investigate fluorescence detection using a standing wave of blue-detuned light of one or more atoms 
held in a deep, microscopic dipole trap. The blue-detuned standing wave realizes a Sisyphus laser cooling 
mechanism so that an atom can scatter many photons while remaining trapped. When imaging more than 
one atom, the blue detuning limits loss due to inelastic light-assisted collisions. Using this standing wave probe 
beam, we demonstrate that we can count from one to the order of 100 atoms in the microtrap with sub-poissonian 
precision. 


PACS numbers: 37.10.De,37.10.Gh,32.50.+d,42.50.Lc 

I. INTRODUCTION 

Atomic ensembles stored in microscopic dipole traps pro¬ 
vide an excellent system in which to study atomic collisions 
and light scattering processes. These microscopic dipole 
traps, so-called microtraps, enable the generation of atomic 
samples of very high density at modest atom number and tem¬ 
perature. Such a system holds the potential for generating mi¬ 
croscopic Bose Einstein condensates [1]. Additionally, the 
use of Rydberg states within these small and dense atomic en¬ 
sembles has applications in quantum information, such as the 
realization of collective qubits [2] and non-linear absorption 
and dispersion of light through collisions between Rydberg 
atoms [3, 4]. 

One of the key challenges in this system is the determina¬ 
tion of the number of atoms in the microtrap. In particular, 
determining the atom number below the level of poissonian 
noise is important for realizing collective quantum gates based 
on Rydberg “superatoms” [2, 5]. A standard approach is to 
fluorescence image the atomic sample by exposing it to red- 
detuned optical molasses and detecting a portion of the scat¬ 
tered light through high numerical aperture (NA) optics [6, 7]. 
The laser cooling provided by the optical molasses enables an 
atom to scatter many photons while remaining trapped. How¬ 
ever, the optical molasses drives light-assisted collisions that 
cause rapid pair-wise loss from the trap; this occurs because 
the red-detuned light excites a colliding pair of atoms onto 
an attractive molecular potential energy curve that is typically 
many times deeper than the optical trap [8]. This technique 
has been used to produce single atom occupancy in a micro¬ 
trap [6, 9], but is not well-suited to the non-destructive count¬ 
ing of more than one trapped atom. One can avoid this loss 
by releasing the atoms from the dipole trap and counting the 
atoms once the density is sufficiently low, either in free-space 
[10], or by recapturing them in a larger volume trap such as a 
magneto-optical trap (MOT) and probing them there [11-14]. 
Recently, atom counting in a MOT has been demonstrated 
with single atom resolution for over 1000 atoms [15]. How¬ 
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ever, if the atomic sample is to be further manipulated within 
the microtrap, or if one wishes to count atoms in several adja¬ 
cent microtraps, it would be particularly useful to employ an 
in situ method. 

Here, we describe an in-trap fluorescence detection method 
that allows us to count one to more than a hundred atoms held 
in a microscopic dipole trap, thereby permitting subsequent 
probing of the trapped atomic cloud. This is an extension 
of earlier work [16], where fluorescence was induced using 
a standing wave of light that was blue-detuned from atomic 
resonance; this configuration realizes a laser cooling mecha¬ 
nism that works well for small detunings ((5 ~ 27r x 15 MHz) 
and for intensities in excess of saturation, thereby yielding 
high photon scattering rates that make it ideal for fluorescence 
imaging. The sign of the detuning and the spatial variation 
of intensity lead to a form of Sisyphus cooling [17, 18] that 
allows an atom to remain trapped during imaging and be de¬ 
tected with near unit efficiency. This capability permits us to 
eliminate the artifacts of loss during the imaging process. For 
samples containing more than one atom, the small blue detun¬ 
ing limits the losses from inelastic light-assisted collisions, 
both by reducing their rate via optical shielding, and, in the 
case where such an inelastic collision occurs, by limiting the 
energy gained by the colliding atom pair to the detuning of the 
imaging light hS [19, 20]. We show that this method can count 
atoms in samples of very high density (n ~ 5 x 10'^ cm“^) 
with sub-poissonian precision. 

The paper is organized as follows. In the next section we 
describe the experimental setup. In section III, we summarize 
the methods used to calculate the light shifts of the dipole trap 
and probe light. Section IV describes fluorescence detection 
of a single atom held in a microtrap, and section V extends 
this to the case of multiple atoms. In section VI we present 
conclusions and provide an outlook for future work. 

II. EXPERIMENTAL SETUP 

We prepare small atomic samples by loading a microscopic 
dipole trap from a cloud of ®®Rb atoms held in a magneto¬ 
optical trap. Figure 1 is a schematic of the setup. The MOT 
operates on the F = 3 to F' = 4 transition of the D2 line 
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FIG. 1. (Color online) Schematic of the experimental setup. ®®Rb 
atoms are loaded from a MOT into a microscopic dipole trap. The 
microtrap is formed by passing 828 nm light through a high NA lens 
mounted inside the vacuum chamber. Atoms held in the microtrap 
are probed by a standing wave of light on the D1 line of ®®Rb. The 
MOT cooling beams are not shown. A portion of the light scattered 
by the trapped atoms is detected by an EMCCD camera. 


of ®®Rb (5 ^Si/ 2 to 5^P3/2) at wavelength A = 780 nm. The 
light for the MOT is supplied as three retro-reflected beams 
with gaussian intensity profiles, each with spot size (1/e^ 
waist) 6.4 mm and apertured to 7 mm in diameter; through¬ 
out the paper, we state the power measured after the apertures. 
A compressed MOT (CMOT) phase typically lasting 150 ms 
overlaps the center of the MOT cloud with the dipole trap. 
Subsequently, the magnetic field is extinguished and 5 ms of 
polarization gradient cooling leads to the loading of approx¬ 
imately one hundred atoms in the microtrap. The number of 
loaded atoms is typically varied by changing the duration of 
the MOT loading time (~ 1 s) and the CMOT phase. Three 
sets of Helmholtz coils serve to null the ambient magnetic 
field. 

The microtrap is formed by focusing a linearly polarized 
gaussian beam with wavelength 828 nm by a high numeri¬ 
cal aperture lens held in the vacuum chamber to a spot size 
of wq = 1.8 pm. The lens is a single-piece molded asphere 
with numerical aperture 0.55 and working distance 2.92 mm 
(Lightpath model 352230). The size and gaussian inten¬ 
sity distribution of the dipole trap beam were confirmed in a 
replica setup outside the vacuum chamber. These dipole beam 
parameters generate microscopic optical traps that are several 
mK deep using ^ 10 mW of optical power; a detailed de¬ 
scription of the light shifts induced by the dipole trap beam is 
presented in section III. 

To detect atoms in the microtrap, we induce fluorescence 
with a probe beam on the D1 line of ®®Rb. Typically, the probe 


beam is blue-detuned from the F = 2 to F' = 3 transition on 
the D1 line (5^Si/2 to 5^Pi/2) for an atom at the center of the 
dipole trap. We choose the D1 line for imaging so that we can 
block stray light from the MOT and dipole trap using standard 
optical filters. The detuning is given by <5 = w — wg, where 
u! is the frequency of the laser field, and ujq the frequency of 
the atomic resonance; a blue detuning corresponds to <5 > 0. 
The probe beam detuning an atom experiences depends on 
its position in the microtrap. In this paper, we reference the 
probe detuning to that experienced by an atom at the center of 
the dipole trap, and denote this by Sc- The probe beam is lin¬ 
early polarized (along direction y in Fig. 1) and retro-reflected 
to form a standing wave by a metal mirror placed outside the 
vacuum chamber. A mechanical shutter lies in the return beam 
path so that the probe beam can be toggled between travel¬ 
ing and standing wave configurations within the course of an 
experimental run. The probe beam has a gaussian intensity 
profile and is focused with a waist of wq = 150 pm at the 
position of the atoms. We use a focused beam primarily to 
avoid it hitting the high NA lens so as to minimize stray light 
in images. Due to the available ports of the vacuum chamber, 
the probe beam runs at an angle of 9° to the radial axis of the 
dipole trap. 

Given there is no cycling transition on the D1 line, we re¬ 
quire repump light to return the atoms to the F = 2 ground 
state. For this, we primarily use the MOT cooling beams de¬ 
tuned by an amount 6^^ from the D2 F = 3 to F' = 3 tran¬ 
sition for an atom at the center of the trap, averaged over the 
magnetic sublevels and assuming a uniform population distri¬ 
bution across these levels. We may also use the F = 3 to 
F' = 2 transition on the D1 line to optically pump the atoms 
back to the lower ground state. This beam shares the same 
optical path as the probe beam; i.e., both D1 beams propagate 
through the same optical fiber. Since our imaging system de¬ 
tects all D1 light, this repump beam increases the signal by 
about 80% [16]. For imaging, the D1 and D2 beams are ap¬ 
plied with rectangular pulse envelopes. 

A portion of the scattered light by the atomic sample is 
collected by the high NA lens and directed onto an Electron 
Multiplying Charge Coupled Device (EMCCD) camera. As¬ 
suming fluorescence into a solid angle of Ttt, the high NA 
lens collects 10% of the scattered light, and half of this is 
reflected by a Polarizing Beam Splitter (PBS). The remain¬ 
ing light is optically filtered for dipole trap light at 828 nm 
and D2 light at 780 nm. Room light is further suppressed 
using shields that enclose the vacuum chamber and surround¬ 
ing optics. An image is formed on the EMCCD by a high 
quality lens with focal length 200 mm. The combined trans¬ 
mission of the PBS, filters and additional optics is 37%. The 
EMCCD has a measured quantum efficiency of rj = 60%, so 
that the total detection efficiency is 2.3%. We typically oper¬ 
ate the camera with multiplication gain set such that the total 
(mean) gain G = yg = 20, where the number of analog-to- 
digital units (ADU) and the incident number of photons are 
linked by A^adu = GN-p^. The gain of the EMCCD is a re¬ 
sult of several amplification stages in series leading to an addi¬ 
tional multiplicative noise term of approximately x/2 [21, 22]. 
The imaging system has magnification 44, and the camera has 
























3 


pixel-size 16x 16 /im^. The image of the atomic sample is in¬ 
tegrated over an 11x11 pixel region of interest, from which 
we subtract the mean background level from two 80 x 80 
pixel) regions on either side of the atomic signal [23]. 

The number of atoms held in the trap can be reduced by en¬ 
gaging the probe beam further blue-detuned from resonance 
[8]. This technique excites an atom pair onto a repulsive 
molecular potential curve; the energy gained by the atom 
pair in the process is at maximum hS, in contrast to the case 
of excitation by a red-detuned beam. In previous work, we 
have achieved a single atom loading efficiency in excess of 
90% [24, 25]. In this paper, we study the imaging of single 
atoms produced by this technique, but we also employ it to 
obtain a single atom fluorescence signal as a reference level 
when imaging multiple trapped atoms. 


III. LIGHT SHIFTS 

The use of a deep dipole trap to hold an atomic sample 
means that the atoms can scatter many photons and remain 
trapped, but it also requires us, when imaging, to take account 
of the significant light shifts induced by the microtrap [26]. 
Furthermore, the blue Sisyphus mechanism used for probing 
relies on the spatial variation of the light shifts induced by the 
standing wave probe beam, so an accurate calculation of these 
shifts is critical to analyzing this mechanism applied to in¬ 
trap fluorescence detection. Our approach is to calculate the 
light shifts generated by the dipole trap and to include these 
as effective spatially dependent detunings in the calculation of 
the light shifts due to the probe beam. This approach is valid 
and convenient because the large detuning from any optical 
transition of the dipole trap light means that the dressed states 
within the dipole trap are very well approximated by the bare 
atomic states with shifted energy levels [27]. 

We calculate the light shifts induced by the dipole trap in 
a semi-classical picture where an applied laser beam is rep¬ 
resented by a classical electric held E(y) cos(a;f). To second 
order in perturbation theory, the shift of atomic level |i) is 
given by [28] 



( 1 ) 

where the sum is taken over all other electronic states |j), and 
Mij (ujij) denotes the dipole matrix element (atomic transition 
frequency) between states |i) and |j). If the frequency of the 
applied held is close to a given atomic transition, the rotating- 
wave approximation may be applied and the second term in 
parenthesis in Eq. (1) may be discarded, leaving the denomi¬ 
nator as the detuning from the chosen transition S = ui — uiij. 

To calculate the light shifts for ®®Rb, we use values for 
the dipole matrix elements and bare atomic transition fre¬ 
quencies [29] to perform the sum in Eq. (1) over the rele¬ 
vant dipole allowed transitions on the D1 and D2 lines [30]. 
The electric held strength is related to the intensity by 
I(r) = iceo \E{v)^. The gaussian intensity distribution of 
the dipole trap beam is /(r) = Jq exp(—2r^/r(;(z)^), where 
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FIG. 2. (Color online) Light shifts for relevant states of ®®Rb along 
the radial direction of the microtrap at a power of 20 mW. (a) States 
addressed by D1 probe beam: Ground state F — 2 (black line) and 
excited state F' — 3 (red line). The probe beam optical frequency u 
has been subtracted from the excited state energy. is the detuning 
of the probe beam at the center of the microtrap. Away from the cen¬ 
ter, an atom experiences a position-dependent detuning S{x) that is 
larger than Sc- (b) States addressed by the D2 MOT cooling beams: 
Ground state F = 3 (black line) and excited states F' = 3 (red 
line) and F' — 4- (blue line). The MOT cooling beams optical fre¬ 
quency uJo 2 has been subtracted from the excited state energies. The 
magnetic sublevels in the excited states on the D2 line experience 
different light shifts; the curves show the mean across the magnetic 
sublevels. Sc^ denotes the detuning of the MOT cooling beams from 
the F = 3 ^ F' = 3 transition on the D2 line at the center of the 
microtrap, with 5^^ (x) the position dependent detuning. 


Iq = 2P/{ttw{z)‘^), P is the beam power, and the position 

1 

dependent spot size is given by w{z) = wq [l -F (^/zr)^] 
with 2 r = ttwq /A the Rayleigh range. 

Eigure 2 shows light shifts calculated using Eq. (1) for rel¬ 
evant states in ®®Rb along the radial direction of the dipole 
trap beam. The red detuning (S < 0) of the dipole trap 
beam shifts the ground states down and the excited states 
up. Eigure 2(a) shows the light-shifted energies of the states 
addressed by the D1 probe beam. Due to the use of lin¬ 
early polarized light, the magnetic sublevels of the ground 
states E = 2,3 and the excited states F' = 2,3 on the 
D1 line are shifted equally. We denote the energy of the 
lower ground state U (r) and the E' = 3 excited state on the 
D1 line ^(r). Eor 20 mW power in the dipole trap beam, 
[14(0) — U{0)]/h = 2 tt X 82.7MHz. The microtrap leads to 
a position dependent detuning i5(r) = <5 — [I4(r) — U{r)]/h 
for the probe beam. Eigure 2(b) shows the light-shifted en- 
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ergies of the states addressed by the MOT cooling beams on 
the D2 line; The F = 3 ground state and the F' = 3,4 ex¬ 
cited states. The light shifts for the excited states on the D2 
line are rriF-level dependent; the curves show the mean over 
the magnetic sublevels, assuming a uniform population dis¬ 
tribution across them. As for the probe beam, the microtrap 
leads to a position dependent detuning for the MOT cooling 
light from the F = 3 to F' = 3 transition, <5^^ (r). However, 
on the D2 line, the situation is complicated by the fact that as 
the atom moves away from the center of the trap, the F = 3 
to F' = 4 transition also becomes near-resonant. 

To calculate the energy shifts of the atom due to the probe 
beam we use the dressed atom approach [17]. For fluores¬ 
cence imaging, the atoms are exposed to a standing wave of 
light that is near-resonant with the F = 2 to F' = 3 transi¬ 
tion on the D1 line. Given the small detuning, we approxi¬ 
mate the system as a two-level atom exposed to a single mode 
light field, described within the electric dipole and rotating- 
wave approximations; i.e., the Jaynes-Cummings model [31]. 
To account for the light shifts induced by the dipole trap, we 
use the position dependent probe detuning S(r) in the calcu¬ 
lation of the dressed states. Our standard imaging conditions 
use (5c ~ 27r X 15 MHz and an intensity higher than the on- 
resonance saturation intensity. The dressed state energies of 
the combined light-atom system are [17] 

^ 2 n(r) = (n++ U(r) - (2) 

where n is the number of photons in the light 
field. The generalized Rabi frequency is given by 

H(r) = [(5(r)2-b a;^(r)] , with the on-resonance Rabi 

frequency wi(r) in the semi-classical picture defined by 
• F(r)//i, where <f>(r) is the phase of the 
light-atom coupling. 

Figure 3 shows the dressed state energies for a ®®Rb atom 
held in the dipole trap, exposed to the standing wave probe 
beam with power 30 fj,W and (5c = 27r x 15 MHz. For sim¬ 
plicity, we treat the probe beam as propagating orthogonal to 
the dipole trap beam. With the choice of a small blue detuning 
for an atom at the center of the trap, the effective detuning of 
the probe beam increases away from the center. Thus, at the 
center of the trap, the effective probe beam detuning is mini¬ 
mal and the light shift induced by the probe beam is maximal. 
Away from the center, the effective blue detuning increases 
and the probe beam light shifts decrease in magnitude. 


IV. SINGLE ATOM DETECTION 

Our main goal in single atom fluorescence detection is to 
detect as much light scattered by the atom as possible while 
leaving it trapped. As such, the two quantities we use to char¬ 
acterize the detection performance are the integrated fluores¬ 
cence - the number of ADUs we detect on the camera - and 
the retention probability - the probability that a single atom 
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FIG. 3. (Color online) Calculated dressed state energy for an atom in 
the microtrap dressed by a standing wave of near-resonant light. The 
probe beam has power 30 fiW and 5c = Stt x 15 MHz. The ground 
state F — 2, mp = 0 is shifted down by the dipole trap (black 
dashed line), and modulated by a standing wave of blue-detuned light 
on the D1 line (black solid line). The excited state F' — 3, m'p — 0 
is shifted up by the dipole trap (red dashed line), and modulated by 
the standing wave (red solid line). 


remains in the microtrap following a fluorescence image. The 
retention probability is a measure of the heating caused by the 
detection process. The higher the single atom scattering rate, 
the higher the integrated fluorescence signal and the greater 
the accuracy with which we can distinguish between the pres¬ 
ence and absence of an atom in the microtrap. The integrated 
fluorescence and the retention probability are coupled in the 
sense that a low retention probability leads to the loss of an 
atom from the microtrap during the detection process, and 
thereby a lower signal than if the atom had remained trapped 
throughout the exposure. 

Figure 4(a) is a schematic of the experimental sequence 
we use to investigate fluorescence detection of a single atom 
held in the microtrap. The sequence begins by loading ^ 
100 atoms into the microtrap from the MOT, whereupon we 
drive inelastic light-assisted collisions to reduce the number of 
atoms in the trap to one or zero [24, 25]. The first fluorescence 
image confirms the presence of an atom. For this image, the 
probe beam is applied as a standing wave with 30 fiW power 
and detuning 6c = 2 tt x 15 MHz, and the D2 beams each has 
power 1.4 mW and detuning 6^^ = 27r x 2.8 MHz. By this 
imaging procedure, an atom is retained in the microtrap with 
99% efficiency. In the second ‘variable’ image, an imaging 
parameter is varied, and the integrated fluorescence for this 
setting is recorded. The third image, identical in parameters 
to the first, tests whether the atom is still trapped after the 
second image. The outcome of the third image, conditioned 
on there being an atom present in the first image, gives the 
retention probability. All images have 10 ms exposure time. 
Unless otherwise stated, each data point shows the mean of 
200 (unconditioned) experimental runs and error bars show 
the standard error of the mean; to avoid obtaining unphysical 
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FIG. 4. (Color online) (a) Experimental sequence to investigate the 
fluorescence detection of a single atom, (b) Fluorescence spectrum 
observed when the probe beam detuning is varied. Integrated flu¬ 
orescence conditioned on the detection of a single atom in the first 
and third image (black circles), data conditioned only on there be¬ 
ing an atom in the first image (red squares). Solid lines are fits of 
asymmetric Lorentzian functions to the data. 


retention probabilities exceeding 100%, we do not correct for 
the small loss induced by the first and third images when cal¬ 
culating the probability of retaining an atom. With a view to 
presenting qualitative arguments to describe single atom de¬ 
tection, we simplify the physical picture by omitting the D1 
repump light on the F = 3 to F' = 2 transition. The addition 
of the D1 repump light does not change the performance of 
the fluorescence detection other than to increase the observed 
signal, a feature we make use of when imaging samples con¬ 
taining several atoms (see section V). 

To investigate the characteristics of the blue Sisyphus 
mechanism as a function of detuning, and to confirm the light 
shift calculations of section III, we performed an in-trap flu¬ 
orescence spectrum of a single atom. To do this, we var¬ 
ied the detuning of the probe beam in the ‘variable image’ 
in the experimental sequence (see Fig. 4(a)), while holding 
the power at 30 /iW. The MOT cooling beams were set to 
= — 27r X 5 MHz and 1.4 mW per beam. Figure 4(b) 
shows the observed single-atom spectrum. The red squares 
show the mean integrated fluorescence for the experimental 
runs where there was an atom present in the first image. The 
black circles show the mean integrated fluorescence condi¬ 
tioned on there being an atom present both in the first and third 
images. The two data sets have been fitted with an asymmetric 
Lorentzian model [32]. The unconditioned data is asymmet¬ 


ric because the standing wave probe beam has a cooling effect 
on the blue side of the spectrum and a heating effect on the 
red side. This sign dependence of the detuning also shifts the 
peak of the unconditioned data approximately 3 MHz to the 
blue. In contrast, the peak of the conditioned data agrees well 
with Sc = 0, confirming the validity of the light shift calcu¬ 
lation. This illustrates a great advantage of non-destructive 
imaging at the single atom level: It allows one to identify ar¬ 
tifacts from features. In this case, we avoid the corruption of 
the spectroscopic signal due to atom loss. 

The Sisyphus mechanism realized by the standing wave 
probe beam cools an atom for blue detunings and heats for red 
detunings [17], so it is somewhat surprising that for modest 
red detunings of the probe beam, —27r x 15 < (5c < 0 MHz, 
an atom remains trapped in some of the experimental runs. 
There are two reasons for this, both relating to the position 
of the atom in the microtrap and the effect this has on the 
frequency of the applied light fields for probing and optical 
pumping. First, the effective probe beam detuning is pushed 
to the blue as an atom travels into the wings of the micro¬ 
trap. As such, although the standing wave probe beam heats 
an atom at the center of the trap for Sc ~ —27r x 15 MHz, in 
the wings the effective blue detuning is restored and the atom 
experiences the blue Sisyphus cooling mechanism, thereby re¬ 
maining trapped. Second, the MOT cooling beams laser cool 
a trapped atom when it moves to the wings of the microtrap: 
Although tuned so as to optically pump an atom at the center 
of the trap back to the F = 2 ground state, i.e., Sc^ « 0, 
away from the center of the microtrap, an atom experiences 
the MOT cooling beam frequency closer to the F = 3 to 
F' = 4 cycling transition, thereby realizing a standard red- 
detuned optical molasses. As such, an atom can be heated at 
the center of the trap by the probe beam, but be cooled by the 
MOT beams in the wings of the microtrap so that it remains 
trapped. 

This shows that it is important to characterize the role the 
MOT cooling light plays in the detection process, even though 
its nominal role is to optically pump a trapped atom to the 
electronic ground state addressed by the probe beam. To do 
this, we varied the detuning of the MOT cooling light while 
keeping the D1 imaging parameters constant. The experiment 
was performed for standing and traveling wave probe beams 
in the second image in order to test the cooling effect of the 
MOT cooling light in the presence and absence of the Sisy¬ 
phus cooling mechanism. As noted above, the blue Sisyphus 
cooling mechanism relies on the spatial modulation of the 
probe beam intensity: In a traveling wave configuration, the 
probe beam has a flat intensity distribution over the extent of 
the microtrap and therefore provides no cooling. Figure 5(a) 
shows the probability of detecting an atom in the third image, 
conditioned on it being detected in the first. For both stand¬ 
ing and traveling wave detection and Sc^ > — 27r x 5 MHz, 
the retention probability is in excess of 97%. However, below 
this detuning, both standing and traveling wave probe beams 
show an approximately linear decrease in retention probabil¬ 
ity, but the decrease is much sharper for the traveling wave 
case. Thus, for S^^ > —27r x 5 MHz, the D2 beams perform 
enough laser cooling to keep an atom trapped while imaging 
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FIG. 5. (Color online) Imaging performance as a function of MOT 
cooling light detuning from F = 3 to F' = 3, (5c ^ , at constant MOT 
cooling light power. The experiment was performed using standing 
wave (black circles) and traveling wave (red squares) probe beams, 
(a) Probability of retaining an atom in the microtrap, (b) Integrated 
fluorescence. 


for this choice of probe beam parameters, whereas below this 
detuning, the blue Sisyphus cooling mechanism of the probe 
beam is necessary to keeping the atom trapped. Figure 5(b) 
shows the integrated fluorescence counts of the second im¬ 
age. In the standing wave case, the counts are maximized for 
' = -27r X 5 MHz. 

To investigate this imaging process further, we scanned the 
probe beam power for several values of detuning. In order 
to focus on the cooling effect of the blue Sisyphus mecha¬ 
nism, the D2 light was set to (5p = —27r x 10 MHz with 
1.4 mW per beam. This value of the detuning was a compro¬ 
mise between accentuating the cooling effect of the standing 
wave probe beam, and ensuring efficient optical pumping and 
hence a good signal. Figure 6 shows the retention probability 
and integrated fluorescence for three values of the probe beam 
detuning: Sc = 2ttx 5, 15 and 35 MHz. 

The retention probability shows an interesting structure for 
the low detuning Sc = 2tt x 5 MHz, which is approximately 
equal to the natural linewidth Fdi = 27r x 5.75 MHz of the 
transition. For zero power, the retention probability is of 
course approximately one, limited only by the performance of 
the first and third images. However, as the power of the probe 
beam is increased very slightly, the retention probability drops 
sharply and reaches a local minimum at a power of 2.5 /rW, 
corresponding approximately to saturation intensity /sat for 
the transition. As the power is further increased, the retention 
probability increases to a peak value and eventually begins to 




FIG. 6. (Color online) Single atom detection as a function of probe 
beam power for three detunings 5c = Svrx 5 (black circles), 15 (red 
squares) and 35 (blue diamonds) MHz. (a) Probability of retaining 
an atom in the microtrap, (b) Integrated fluorescence. 


decrease. The initial drop at low power can be understood 
as Doppler heating predominating over the Sisyphus mecha¬ 
nism. To give a qualitative picture of this behaviour, we argue 
using the simple treatment of Doppler heating/cooling, where 
the total force on an atom due to a counter-propagating pair of 
laser beams is given by the sum of the contribution from each 
beam [33]. Using this description, the Doppler heating rate 
attains its maximum value for Sc = Fdi and / = /sat. Above 
2.5 pW, the Doppler heating rate decreases and the cooling 
effect of the blue Sisyphus mechanism increases as the mag¬ 
nitudes of the dressed state energies increase at the anti-nodes 
of the standing wave. The minimum temperature achievable 
with the blue Sisyphus mechanism is proportional to the mod¬ 
ulation amplitude of the dressed states E^od in the standing 
wave [34]; as the power is further increased, the temperature 
of the sample increases and eventually becomes commensu¬ 
rate with the trap depth, so that the D2 cooling beams can no 
longer provide enough cooling to keep the atom trapped. If 
one assumes a Maxwell-Boltzmann distribution for the energy 
of an atom with temperature T ^ E^od / Kb in the microtrap, 
integrating the distribution from zero to Uq shows the same 
trend as this experimental retention probability. 

At higher blue detunings, the dip in retention probability 
at low power is significantly reduced: Sc = 27t x 15 MHz 
shows a small decrease at 10 ^W, and at Sc = 2'k x 35 MHz 
the retention probability is essentially flat for low power. This 
can be understood from the fact that for detunings larger than 
the natural linewidth, the peak Doppler heating rate shifts to 
higher intensities, where the modulation of the dressed state 
energies (Eq. (2)) is already significant, and thus the cool¬ 
ing effect of the blue Sisyphus mechanism predominates over 
Doppler heating. 
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V. MULTI-ATOM DETECTION 


In Fig. 4, we saw that trap loss can distort the fluorescence 
signal of a single atom. However, due to the fact that one 
atom can be detected with near-unit efficiency, this loss effect 
could be eliminated by adding additional pictures to the ex¬ 
perimental sequence and appropriate post-processing. In the 
case of fluorescence detection of multiple atoms in a micro¬ 
scopic volume, the process is complicated by rapid two-body 
loss due to light-assisted collisions between atoms. In previ¬ 
ous work [16], we showed that fluorescence imaging using a 
blue-detuned standing wave of light could be used to count 
three atoms held in a microtrap based on the unambiguous 
signature of steps in the integrated fluorescence. However, 
this signature is lost when counting larger atomic samples be¬ 
cause light-assisted collisions reduce the atom number during 
a given exposure time, thereby blurring the steps in fluores¬ 
cence. Thus, to accurately count multiple atoms in a micro¬ 
scopic volume, we must devise additional techniques to miti¬ 
gate the effects of loss during the imaging sequence. 

In our system, inelastic light-assisted collisions during flu¬ 
orescence imaging are the dominant loss mechanism. In gen¬ 
eral, binary collisions between cold alkali atoms in their elec¬ 
tronic ground states occur through an induced dipole-dipole 
interaction, a second order effect with interaction energy given 
by Cq/R^, where R is the inter-atom distance and Cq is the 
van der Waals interaction coefficient [ 8 ]. In the presence of 
near-resonant light, the electronic excited state population is 
non-zero and atoms can interact via a resonant dipole inter¬ 
action Cs/R^. One can estimate the inter-atom distance at 
which the light is resonant with the two-atom system, the 
Condon point Rq, by hS = Cs/Rq. A typical C 3 value is 
20 eVA3, so for b = 27r x 15 MHz, Rc - 500 A. For 
100 atoms in the microtrap at 160 pK, the peak density is 
no ~ 5 X 10^^ cm“^, and the mean inter-atom distance is 
R ^ -250 A. Thus, in this system the mean inter¬ 

atom distance is of the same order as the interaction length, 
and light-assisted collisions predominate. 

Our approach to counting atoms in mesoscopic samples is 
to model the effect of two- and one-body loss on the time- 
resolved integrated fluorescence and from this determine the 
atom number as a function of imaging time. A similar ap¬ 
proach using standard red-detuned cooling light was recently 
demonstrated for up to — 20 atoms in a larger volume dipole 
trap [5]. We model the atom loss during imaging with a mas¬ 
ter equation, in which the time evolution of the probability p n 
to obtain N atoms in the microtrap is given by [35, 36] 


dpN 

dt 


/3(E2 -1) 


N{N- 1) 


-PN 


-7(E-l)[7Vp^], (3) 


where [3 and 7 are the two-body and one-body decay rates, re¬ 
spectively, and E is the step operator. This operator is defined 
by its effect on a function f{N) by E[/(A^)] = f{N + 1) and 
E~^[/(iV)] = f{N — 1). The first term in Eq. (3) describes 
a two-body loss process, where both atoms are lost from the 
trap; we call this “process 2”. One can also model a two-body 
collision where only one atom is lost from the trap (“process 
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FIG. 7. (Color online) Atom number distributions evaluated for pro¬ 
cesses 1 and 2 from Eq. (3) using Nq = 63, /3 = 4.9 s“^ and 
7 = 0.97 s~^. To ensure the same mean atom number for the two 
processes, we use /3 for process 2 and 2/3 for process 1 (see Eq. (4) 
and associated text), (a) Mean atom number (black line) calculated 
from Eq. (4) assuming /3AA^ = 0. The shaded regions indicate 
the 2(t width of the atom number distributions for process 2 (gray) 
and 1 (red), (b) Atom number distribution at 16 ms for process 2. 
The solid line shows the poissonian distribution for the mean atom 
number N = 10.6. (c) Same as (b) but for process 1. 


1” ) by replacing E^ in this term by E; in our experiments, 
such a process arises from the combined effect of the imag¬ 
ing beam and the MOT cooling beams used to optically pump 
the atoms. By probing an atomic sample containing just two 
atoms, we And that process 1 and process 2 occur with ap¬ 
proximately equal probability. Equation (3) comprises a sys¬ 
tem of equations that may be solved numerically for samples 
containing several hundred atoms. 

The time evolution of the mean number of atoms, 
N = X]w=o d^PN, may be found from Eq. (3), yielding 

^ = -j3N(N-l)--fN-l3AN^, (4) 

dt 

where AN'^ is the variance in the atom number. Equation (4) 
also describes process 1 , but in this case the first term on the 
right hand side reads —^N{N — 1 ), given that only one of the 
two atoms is lost from the trap. Accordingly, to achieve the 
same mean decay rate when comparing the two processes, we 
use decay rate /3 for process 2, and 2/3 for process 1. 

Eigure 7 shows the time evolution of the atom number dis¬ 
tribution calculated using Eq. (3) for representative parame¬ 
ters used in our experiments. Beginning from 63 atoms in 
the trap, the imaging leads to a decrease in the mean num¬ 
ber of atoms and an increase in the atom number variance. 
Eigure 7(a) shows TV as a function of exposure time and the 
shaded regions denote the 2 cr width of the atom number dis- 
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tribution for processes 1 and 2. The solid line shows the 
mean atom number calculated from Eq. (4) and assuming 
/3AiV^ = 0; it is clear that this solution approximates the 
means of the two distributions very well, showing that this 
correction to the mean atom number is negligible for our ex¬ 
perimental parameters. The atom number variance initially 
grows from its (idealized) value of zero at f = 0 to reach its 
“limiting” value at f « 10 ms. This value is determined by 
the competition between the one and two-body loss mecha¬ 
nisms: Acting alone, the former leads to a poissonian number 
fluctuations, i.e., AN"^ = N, whereas process 2 induces an 
atom number variance of 0.667iV for IV > 10 [37], reducing 
to ATV^ = 0.57V for the special case of TV = 0.5 [6]. Process 
1 , on the other hand, leads to zero atom number variance in 
the limit that TV —1, because the reduction in atom number 
stops when there is only one atom left in the trap. The dif¬ 
ferences in the atom number variance produced by these two- 
body processes are evident in Figs. 7(b) and (c), which show 
the atom number distributions at f = 16 ms for process 2 and 
1, respectively. Process 1 induces a smooth, unimodal distri¬ 
bution, whereas process 2 yields a jagged distribution with a 
unimodal envelope. For this calculation, odd atom numbers 
are more probable because the initial atom number was cho¬ 
sen to be odd, and process 2 leads to pairwise loss. In both 
cases, however, the width of the distribution is smaller than a 
poisson distribution with the same mean. 
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FIG. 8 . (Color online) (a) Experimental timeline for multi-atom 
detection, (b) Integrated fluorescence as a function of exposure 
time for three values of the initial atom number TVq. Each data set 
was fitted with Eq. (7) using Fi and 7 as fixed input parameters 
(see text). The fitted parameters were: (black circles) Ao=85.9(5), 
f} = 6.41(4)s“^; (red squares) iVo=58.4(4), /3 = 6.00(5)s“^; (blue 
diamonds) TVo=31.3(3),/3 = 5.70(9)s“^. The number in parenthesis 
gives the uncertainty from the fitting routine in the least significant 
figure. 


A. Mean values 


This discussion shows that to determine the mean number 
of atoms in the trap, we may solve Eq. (4) under the condition 
that the atom number variance term is neglected. We will re¬ 
turn to the fluctuations about the mean atom number in section 
V B, where we quantify the precision of the imaging method. 
As noted above, the choice of process 1 or 2 does not affect 
the description of the mean atom number other than the fitted 
value of p. To simplify the comparison of our method with 
standard fluorescence detection methods, we model the mean 
value data with process 2. The solution of Eq. (4) within the 
approximation that /3ATV^ = 0 is given by 


N{t) 


^At 


4(e^*-l) + 


1 

No 


(5) 


where A = /3 — 7 and TVq is the initial atom number. In 
an experiment, we detect the amount of light scattered by the 
atomic sample 


^int(T) 


FiN{t')dt\ 


( 6 ) 


where we have assumed that the integrated fluorescence is 
given by the product of the single atom fluorescence rate Fi, 
which is assumed to be constant in time, and the number of 
trapped atoms. Performing the integration, we obtain 



r/3TVo(e^* 

A 



(7) 


Thus, to determine the initial atom number of a trapped sam¬ 
ple, we take in-trap fluorescence images for several values of 
exposure time and use Eq. (7) to model how the integrated 
fluorescence changes in time. 

Figure 8(a) is a schematic of the experimental sequence 
used to determine the number of atoms held in the microtrap. 
We achieve a given initial atom number TVq by varying pa¬ 
rameters in the MOT and CMOT phases; through this loading 
procedure and evaporation, the atomic clouds thermalize to 
a temperature T ~ 160 pK [38]. The time dependence of 
Tint is probed by exposing the trapped sample to several short 
pulses of imaging light and recording the fluorescence on the 
EMCCD camera. For this fast imaging sequence, the camera 
is operated in kinetics mode and the camera chip masked such 
that we can take ~ 20 frames for our chosen region of inter¬ 
est. Each imaging pulse has 1 ms duration, with a 22 ms delay 
between frames. Tint is generated by taking the cumulative 
sum of the integrated signal for each frame [39]. In order to 
measure Fi, at the end of the fast imaging sequence we drive 
light-assisted collisions with the probe beam set to reduce the 
trapped atom number to zero or one [24, 25], whereupon we 
take a 10 ms fluorescence image (“Image 1” in Fig. 8(a)) using 
the same parameters as for the fast imaging sequence. Finally, 
a standard image (“Image 2”) of duration 10 ms determines 
with near unit efficiency whether there was zero or one atom 
held in the microtrap in “Image 1”. Using the final two im¬ 
ages, and averaging over several experimental runs, we obtain 
the fluorescence rate for one or zero atoms, allowing us to 
calculate the background subtracted single atom fluorescence 
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rate Fi and the background subtracted integrated fluorescence 
The use of collision beam parameters to produce 50% 
zeros and ones minimizes the statistical error in the measured 
background subtracted single atom fluorescence rate, mean¬ 
ing that one could also employ the standard approach of us¬ 
ing light-assisted collisions driven by red-detuned molasses 
for this step. 

Figure 8 (b) shows the integrated fluorescence as a func¬ 
tion of exposure time for three values of the initial 
atom number Nq. The experiment was performed for a 
dipole trap power of 30 mW, yielding a trap depth of 
Uq = h X 87.1 MHz = TTb x 4.18 mK, with standing wave 
probe beam parameters (5c = 27r x 15 MHz and power 
40 /iW for both the F’ = 2 ^ F' = 3 and 
F = 3—>^F' = 2D1 transitions, and MOT cooling beam pa¬ 
rameters (5p = — 27r X 6 MHz at 0.6 mW per beam. The 
data has been fitted with Eq. (7), using Fi and 7 as fixed 
parameters. The single atom decay rate was measured to be 
7 = 0.97 s“^ in a separate experiment where a single atom 
was prepared in the microtrap and exposed to the same imag¬ 
ing conditions as we use for iV-atom samples. The initial 
slope of the integrated fluorescence for each data set is given 
by NqFi . For longer exposure times, the fluorescence rate de¬ 
creases as a result of the loss of atoms from the trap during 
imaging. It is evident that Eq. (7) models the experimental 
data very well, yielding measured initial atom numbers of 31, 
58 and 86 for the three data sets. We have verified the method 
for initial atom numbers up to 130, which is approximately 
the highest number of atoms we obtain in the microtrap using 
the described loading procedure. 

The measured values of two-body loss rate are significantly 
smaller than those reported for comparable experiments that 
use red-detuned imaging light. Our low value of /3 arises 
from the fact that the probe beam is blue-detuned with re¬ 
spect to the primary transition used for imaging: The blue 
detuning excites the colliding atoms onto a repulsive molec¬ 
ular potential curve [19], which limits trap loss due to in¬ 
elastic light-assisted collisions by ensuring that the maximum 
energy the atom pair gains is M; for our experimental pa¬ 
rameters, this gained energy is significantly smaller than the 
trap depth. In order to compare the measured values of /3 
with those from other experimental configurations, we cal¬ 
culate the normalized two-body loss rate /3norm = 132'J2V, 
where V = {2ttKbJ '/is the volume occupied by 
the sample at temperature T, ui is the geometric mean of the 
trapping frequencies, and m is the atomic mass [40]. From 
Fig. 8 , we measure a typical value /3 = 6 s“^, which in nor¬ 
malized form reads /3norm = 1-4 X 10“^^ cm^s“^. For Com¬ 
parable experiments using red-detuned fluorescence imaging 
light, the normalized two-body decay rate is typically two or¬ 
ders of magnitude larger [41]. In particular, we And /?noim = 
1.4 X 10“® cm^s“^ for the experiment described in Ref. [5]. 
It is this feature of our imaging technique that permits the de¬ 
termination of atom number in samples of very high density. 

In our treatment, we assume that the single atom fluores¬ 
cence rate Fi does not depend on the atom number. This as¬ 
sumption may lead to systematic errors for very dense sam¬ 
ples. In particular, light scattering may be suppressed due 


to dipole-dipole interactions [42] or if the temperature of the 
atoms is number dependent. The latter would lead to a num¬ 
ber dependent broadening of the transition due to the trap light 
shifts. An atom number dependent temperature could arise 
due to heating from inelastic collisions that release a low en¬ 
ergy such that the colliding atoms are not lost. Inclusion of 
3-body processes in the model (Eqs. (4) to (7)) did not im¬ 
prove the fit in Fig. 8(b), but for higher densities they may 
also play a role. Reference [5] demonstrated that the atom 
number in small samples can be determined from the Rabi fre¬ 
quency of collective Rydberg excitations, which may provide 
a method for detecting potential systematic errors in our pa¬ 
rameter regime. However, in the following section we present 
a noise model for the imaging process which builds on the fit¬ 
ted values of iVo and /3. This model contains no additional free 
parameters and shows very good agreement with the experi¬ 
mental data, thereby supporting the accuracy of the method. 


B. Fluctuations 


In this section, we develop a noise model to character¬ 
ize the imaging method described above, thereby extending 
the method from the mean atom number to the atom number 
distribution within the microtrap. Similar to previous work 
[15, 43], we characterize the noise properties of the imag¬ 
ing method by calculating the two-sample variance from the 
atom number measured by two consecutive imaging pulses: 
cr^ = ivar(iV 2 — -^i)- Contrary to previous work, however, 
we calculate the noise model based on the full solution of the 
master equation, and use no free parameters to fit the noise 
model to the data. We also compare the atom number at the 
end of the first imaging pulse with the atom number at the 
beginning of the second imaging pulse N^. Such an imaging 
procedure is experimentally relevant, allowing one to deter¬ 
mine the atom number in the microtrap, perform an experi¬ 
ment in which the atom number changes, and finally measure 
the number of atoms remaining in the trap. In light of the dis¬ 
cussion of Eq. (3), one may also view the imaging process as 
a method to prepare an atomic sample with an atom number 
defined better than the poissonian limit. 

The experimental sequence used to characterize the noise 
properties of the imaging method is identical to that shown 
Fig. 8(a), except that in the fast imaging sequence we take 
only two images, each of duration r. The atom number may 
be determined from a single imaging pulse by combining 
Eqs. (5) and (7). At the end of the 5th imaging pulse, the 
atom number is given by 



2^ _ g —/3t^inl/F ■ 

1 - ’ 


(8) 


and at the beginning of the pulse by 


= 


A 


g/3finl/F _ 2^ 


jAr _ I 


(9) 


The data set for this section comprises 100 experimental runs 
for each value of r, with the exposure time sampled between 
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FIG. 9. (Color online) (a) Measured atom number at the end of pulse 
1 and 2 for three imaging pulse durations: r = 0.1 ms (black circles), 
T = 1 ms (red squares), r = 3 ms (blue diamonds). The dashed line 
shows N 2 = A^i. (b) Same as (a), but we compare the atom number 
at the end of the first pulse A^i and the number at the beginning at the 
second pulse N 2 inferred from the second pulse. 


0.1 and 3 ms. The value of /3 was determined using the method 
of the previous section (see Fig. 8(b) and associated text), 
where the integrated fluorescence for each setting of r was 
found by averaging the signal over the 100 experimental runs; 
the obtained value was /3 = 4.94 s“^. The initial atom number 
for the data set was gaussian distributed with mean value 62.5 
and a standard deviation of 8.5 due to fluctuations in the load¬ 
ing process; this was obtained from averaging obtained 
from Eq. (9) over the data set. 

The noise properties of the measured atom number changes 
as a function of the exposure time. Figure 9(a) shows the mea¬ 
sured atom number N 2 v Ni for three values of the pulse dura¬ 
tion. For the shortest pulse duration, the signal shows a large 
degree of scatter; in this regime, the shot noise of the detected 
light is the dominant noise source. For longer pulse durations, 
the scatter decreases markedly but the mean atom loss during 
the second imaging pulse causes the data to deviate from the 
line N 2 = iVi. Figure 9(b) shows N 2 v Ni. Given that Ni 
and N 2 are measurements of the same atom number, the data 
lies on the line N 2 = Ni. However, the inferred atom number 
at the start of the second pulse shows a larger scatter than Ni. 
To understand this behaviour, we model the noise in the mea¬ 
sured atom number arising from the shot noise in the detected 
light and the loss of atoms during the imaging process. 

The light shot noise contribution to the atom number 
variance may be found through standard error propagation. 
In particular, the variance in the atom number at the end 
of an imaging pulse due to light shot noise is given by 
= {dNi/dFiat)'^ajp.^^ + {dNi/dFiYa\,^, assuming negli¬ 
gible uncertainty in the values of (3 and 7. A comparison 
of the two remaining terms in the error propagation formula 
shows that the term proportional to is also negligible. 
Thus, {dNi/dFint )"^gives the variance of the atom num¬ 
ber arising from the light shot noise in pulse i. In princi¬ 
ple, = g^r]Npii = pFint (see discussion of the camera 
in Sec. II). However, due to the spatial variation in the pho¬ 
ton scattering rate for an atom moving in the microtrap and 
imaging light, and the excess noise induced by the EMCCD 
camera, = aFim, where a ~ 110 and is determined 
from the histogram of ADU counts for a single atom exposed 
to the same imaging conditions. Despite the light shielding 


of the apparatus, the background signal is also finite so that 
it makes a contribution cTg to the variance of the integrated 
signal: ap^^ = aFint + o-q, where both terms depend on the 
exposure time. 

To model atom loss during the imaging pulses, we per¬ 
form a Monte-Carlo simulation based on the master equation 
of Eq. (3). The atom number distribution at the end of the 
first imaging pulse is obtained by solving Eq. (3) using the 
initial conditions obtained from the experiment. At time r, 
the atom number is sampled and rounded to the nearest inte¬ 
ger to give Ni', this atom number is then used as the initial 
condition for the atom number distribution’s evolution during 
the second imaging pulse. The resulting probability distribu¬ 
tion is sampled to obtain N 2 . This procedure is repeated 1000 
times for each value of r, from which we calculate the two- 
sample variance due to atom loss, (Tloss- To evaluate the noise 
properties of the method when comparing A^i and N 2 , the 
sampled value of N 2 is inserted into Eq. (5) and rearranged 
to yield N 2 , whereupon we evaluate the two-sample variance 

ivar(Al0 - Ni). 

Figure 10(a) shows the two-sample variance as a func¬ 
tion of exposure time for A^i and N 2 . Also shown 
is the two-sample variance predicted by the noise model 
cr^ = -F + CTlosj. As in Fig. 9(a), the scatter of the 
experimental points initially decreases with the exposure time 
and reaches its minimum at r Ri 0.4 ms. The noise model 
shows that the noise contributions from the light and the atom 
loss are approximately equal at this point. For larger t, the 
atom loss becomes the dominant noise source and the two- 
sample variance increases. In a small interval of t values 
around 0.4 ms, the two-sample variance in detected atom 
number lies below the level of poissonian fluctuations. For 
exposure times longer than 3 ms, the two-sample variance 
again decreases due to two effects arising from the evolution 
of the atom number distribution: The first stems from the de¬ 
creased loss rate at late times, meaning that the relative differ¬ 
ence in mean atom number N 2 — Ni is smaller than for low 
T (see Fig. 7(a)); the second arises from the reduced variance 
in A^i and N 2 at late times due to light-assisted collisions (see 
Figs. 7(b) and (c)). 

The two-sample variance associated with A^i and A’® is 
shown in Fig. 10(b). The variance predicted by the noise 
model is given by cr^ = + ^loss- In the absence of 

the atom number variance associated with mean atom loss, the 
light shot noise is the dominant noise source until r r; 1 ms, 
and remains finite at higher t, consistent with Fig. 9(b) where 
it is apparent that the operation of Eq. (9) to obtain N 2 leads to 
increased variance. Nonetheless, for all exposure times above 
T = 0.1 ms, the imaging technique leads to a two-sample vari¬ 
ance well-below the level of poissonian fluctuations. Indeed, 
for 0.75 < T < 3 ms, the two-sample variance lies a factor of 
three to four below the poissonian limit, which is remarkable 
given the very high density of the atomic sample. 

In calculating the two-sample variance for 7Vi and N 2 , the 
use of process 1 in the noise model yields a significantly lower 
value of noise than the case of process 2, whereas for Ni and 
N 2 , the two loss processes give rise to approximately the same 
noise level. This arises from the fact that the two-sample vari- 
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VI. CONCLUSION AND OUTLOOK 
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FIG. 10. (Color online) (a) Two-sample variance of Ni and N 2 
as a function of exposure time. (Black circles) Experimental data, 
error bars show 2cr statistical error, (red dash-dot line) calculated 
light noise (blue dotted line) calculated noise arising 

from atom loss ctloss based on Eq. (3) using process 1, (black solid 
line) calculated total noise + 0 %^ + ctloss using process 

1 , (green dashed line) calculated total noise model using process 2. 
The shaded region indicates the interval where the atom number vari¬ 
ance lies below the level of poissonian fluctuations, (b) Same as (a), 
but for Ni and Nu¬ 


ance for iVi and N 2 is dominated by the mean atom loss, so 
that the contribution from the widths of the atom number dis¬ 
tributions is small. In contrast, by comparing iVi and N 2 we 
remove the contribution of mean atom loss to the two-sample 
variance, whereby the choice of loss mechanism becomes the 
determining factor. As shown in Fig. 7, process 2 leads to 
a broader number distribution than process 1. Thus, process 
2 leads to a larger uncertainty in atom number N 2 , which is 
amplified by the inversion of Eq. (9) to obtain N^. 

The experimental data in Fig. 10(b) is better matched by 
the noise model assuming process 1, implying that this pro¬ 
cess is the dominant loss mechanism. As noted above, how¬ 
ever, an experiment with just two atoms showed that process 
1 and 2 occur with approximately equal probability. A par¬ 
tial explanation for this discrepancy could lie in the heating 
of the sample due to the high rate of inelastic collisions when 
probing many atoms in a dense sample; In [44] we found that 
an increase in sample temperature led to a higher probability 
for process 1 in inelastic collisions that release an energy less 
than the trap depth. Finally, at large exposure times, the exper¬ 
imental two-sample variance lies slightly below that predicted 
by process 1, indicating that there may be effects that are not 
captured by the noise model. 


We have demonstrated a method of in-trap fluorescence 
imaging based on a standing wave of light that is blue-detuned 
from atomic resonance. In general, imaging in red-detuned 
dipole traps is complicated by the multitude of atomic energy 
levels and their different light shifts induced by the dipole trap 
beam(s). Here, given the use of probe light tuned to the blue of 
the uppermost hyperflne level, F' = 3 on the D1 line of ®®Rb, 
the probe has the greatest effect at the center of the trap, mean¬ 
ing that a two-level atom picture is sufficient. By combining 
the standing wave probe beam with a standard set of MOT 
cooling beams, we achieve a single atom detection probabil¬ 
ity of approximately 99% for a wide range of experimental 
parameters. 

We have explored several features of single-atom detection 
using this method including the role of detuning for the probe 
and MOT cooling beams, and the intensity dependence of the 
Sisyphus cooling mechanism. The use of single atoms to char¬ 
acterize the method was important to avoid spurious results 
arising from the single-atom loss caused by the detection pro¬ 
cess. 

The method was extended to the detection of up to ^ 
100 atoms held in the microtrap, at densities exceeding 
10^^ cm“^. When imaging such high-density samples in a 
trap of finite depth, inelastic two-body light-assisted collisions 
comprise the dominant loss mechanism. We mitigated this 
loss by using a blue-detuned probe beam. The effect of the 
small blue detuning is two-fold: It reduces the rate of inelas¬ 
tic light-assisted collisions via optical shielding, and, when 
an inelastic collision occurs, it limits the energy gained by 
the colliding atom pair to the detuning of the imaging light. 
By choosing the trap depth to be much larger than the en¬ 
ergy gained in an inelastic collision, we achieved very small 
two-body loss rates /3 for such a dense sample. Indeed, the ob¬ 
tained values of the two-body loss rate, normalized to take ac¬ 
count of the microscopic volume of the atomic sample, were 
two orders of magnitude smaller than those obtained using 
the standard approach of fluorescence detection based on red- 
detuned molasses light. 

We modelled the effects of one- and two-body loss on the 
sample using a master equation to determine the atom number 
as a function of imaging time. For two-body loss, we consid¬ 
ered two loss “channels”: Process 1, where one atom is lost 
from the trap due to a binary light-assisted collision, and pro¬ 
cess 2, where both atoms are lost. These two processes lead 
to different atom number distributions in the trap, the signa¬ 
ture of which was evident in the noise analysis of the imaging 
method. The method leads to the preparation and determi¬ 
nation of the atom number up to a factor of four below the 
poissonian limit. When combined with internal state prepara¬ 
tion, this method can be used to generate mesoscopic atomic 
samples useful for quantum information processing. 
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